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Dimensional crossover in Sr 2 RuOi within slave-boson mean-field 
theory 
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Abstract. - Motivated by the anomalous temperature dependence of the c-axis resistivity of 
S^RuCU, the dimensional crossover from a network of perpendicular one-dimensional chains to a 
two-dimensional system due to a weak hybridization between the perpendicular chains is studied. 
The corresponding two-orbital Hubbard model is treated within a slave-boson mean-field theory 
(SBMFT) to take correlation effects into account such as the spin-charge separation on the one- 
dimensional chains. Using an RPA-like formulation for the Green's function of collective spinon- 
holon excitations the emergence of quasiparticles at low-temperatures is examined. The results 
are used to discuss the evolution of the spectral density and the c-axis transport within a tunneling 
approach. For the latter a regime change between low- and high-temperature regime is found in 
qualitative accordance with experimental data. 
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Introduction. — Physical effects due to reduced di- 
mensionality are most prominently visible for strongly cor- 
related electrons. While under generic conditions low- 
energy electronic excitations retain their quasiparticlc 
character in three and most likely two dimensions, corre- 
lated electrons in one dimension fractionalize in separate 
collective charge and spin excitations. All solids are three- 
dimensional and low-dimensionality appears through spe- 
cific highly anisotropic electronic structure. In many 
cases, the effective dimensionality of such systems de- 
pends on the difference of energy scales which can lead 
to dimensional crossovers as systems are cooled or sys- 
tem parameters are changed, for example, by applying 
pressure. The change of effective dimensionality is of- 
ten marked through the onset of ordered states which 
are suppressed in low-dimensional systems due to thermal 
or quantum fluctuations. Such changes are well known 
from organic compounds, such as (TMTSF^X which un- 
dergoes a sequence of transitions as function of temper- 
ature, pressure and chemical composition (X). Another 
remarkable feature is the modification of transport prop- 
erties which can depend on the nature of the charge car- 
riers or on a change from coherent to incoherent trans- 
port. An example of such a system is Sr 2 RuC>4 which 
is a quasi-two-dimensional strongly correlated Fermi liq- 
uid best known for its unconventional superconducting 



phase [1]. Sr2RuC>4 has a single-layered perovskite struc- 
ture with tetragonal symmetry. Down to lowest normal 
state temperatures it displays a strongly anisotropic elec- 
trical resistivity: p c / p a b ~ 10' J at T = IK. Although this 
system shows Fermi liquid behavior (p = p + AT 2 ) in all 
directions at low temperatures (T < AOK) , there is an 
extraordinary difference in the temperature dependence 
of the resistivity at higher temperature for inplane and o 
axis transport. While the inplane resistivity is monotoni- 
cally increasing, along the c-axis it shows a regime change 
from low-temperature metallic to high-temperature insu- 
lating behavior around T* « 130/\ [2]. This change is 
not associated with any structural transition or ordering 
phenomenon and has been attributed to a so far unspeci- 
fied correlation effect of the electrons. For the electrons to 
tunnel between layers, sufficient quasiparticle weight and 
lifetime in the planes, respectively, are required. A rapid 
loss of quasiparticle weight when increasing the temper- 
ature was observed in ARPES measurements [3] and at- 
tributed to scattering processes. Motivated by this be- 
havior we propose here a mechanism based on a dimen- 
sional crossover within the basal plane of Sr2Ru04 which 
affects the inter layer transport. Alternative explanations 
based on electron-phonon coupling [4,5] and charge fluctu- 
ations caused by the tunneling electrons [6] have been dis- 
cussed in the literature. The description of the electronic 
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Fig. 1: a) The two-dimensional nature of the Fermi surfaces 
for the d zx and d yz orbitals is lost for temperatures above the 
weak inter-orbital hybridization so that the electrons behave 
like being in one-dimensional bands, b) The hopping structure 
for the d zx and d yz electrons (here shown for a d zx placed in 
the middle) : in lowest order via an O 2p z to a nearest neighbor 
of the same band or as a next order term via 2 O 2p z to a next- 
nearest neighbor of the other band. 



properties of Sr2Ru04 requires three bands crossing the 
Fermi energy with an average filling of 4/3. Their shape 
is well accounted for by considering the 7r-hybridization of 
the 0-2p-orbitals with 4<i-i2g-orbitals of Ru in each ruthe- 
nium oxide layer [7,8]. The d^-orbital lies in the RuC>2 
plane and forms the genuinely two-dimensional 7-band. 
The other two orbitals d yz and d zx forming the so-called 
a- and /3-band do not hybridize with the 7-band as they 
have different parity with respect to reflection at the basal 
plane. In contrast to the 7-band which has practically no 
interlayer coupling, these other two orbitals are also re- 
sponsible for the weak interlayer hybridization [9]. Thus, 
in order to understand the c-axis transport we will focus 
our attention to the d yz - and (^.-orbitals which have a pe- 
culiar hybridization topology within the basal plane. Both 
orbitals form through nearest-neighbor intra-orbital n- 
hybridization an essentially one-dimensional band disper- 
sion along the y-direction (^-direction) for the d. yz -orbital 



(^zz-orbital). The inter-orbital coupling occurs through 
considerably weaker next-nearest-neighbor hybridization 
between d yz and d zx and leads to two-dimensional bands 
which nevertheless keep strong one-dimensional features 
(see fig. la)) with pronounced nesting properties reflected 
in the enhanced incommensurate magnetic spin fluctua- 
tions [10]. We show here that the structure of these or- 
bitals could lead to a crossover behavior between the one- 
and two-dimensional regime at the energy scale associated 
with the next-nearest-neighbor inter-orbital coupling sim- 
ilar to other systems discussed in literature [11]. For this 
purpose we develop a slave-boson scheme on the model of 
the d y2 -(i Z2 ;-orbitals in order to describe the fractionaliza- 
tion of the quasiparticlcs in the one-dimensional system. 
We will show that this method allows us to qualitatively 
account for the emergence of quasiparticle states in the 
hybridized bands in a rather simple manner. This method 
will also be used to give a qualitative discussion of the 
regime change in the interlayer transport which we con- 
sider in the tunneling limit. 



Formalism. — We start with the two independent, 
one-dimensional bands (henceforth labeled v = 1, 2 corre- 
sponding to the x- and y-directions of the d zx - and d yz - 
orbitals, respectively) described by the Hubbard Hamilto- 
nian 
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,t 



(1) 



with t the hopping integral, U the on-site Coulomb inter- 
action and {i,j) u denoting nearest neighbors in the corre- 
sponding direction. The operator c} via creates an electron 
in the orbital v at site i with spin a and n via = cl, ia c vicr . 
Since such correlated electron systems are difficult to han- 
dle we seek an approximation by assuming that Coulomb 
repulsion is sufficiently strong as to suppress essentially 
double occupancy. In this case we may use a slave-boson 
technique which had been developed in the context of the 
t- J- model. By writing the electron operators c V i a in terms 
of new fermionic operators representing the spin degrees 
of freedom (spinons) and bosonic operators for the charge 
degrees of freedom (holons), c uia = f vi Jy vi and stipulating 
that each site is either occupied by a spinon or a holon, 
double-occupancy is strictly inhibited [12]. From this re- 
placement the so-called slave-boson Hamiltonian results, 
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with the Lagrange multipliers \ V i enforcing the constraints 
of one spinon or holon per site. To analyze this Hamilto- 
nian we resort to the following standard approximations. 
We decouple the first term by introducing mean-fields for 
the coherent motion of holons and spinons as 



xfH^O and X 
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respectively, with i, j being nearest-neighbor indices. The 
treatment of the Hamiltonian is then further simplified by 
taking the mean-fields as being independent of band and 
site and replacing the local constraints by a global one. 

With this approximation we arrive at a slave-boson 
mean-field Hamiltonian representing two independent 
species of particles, spinons and holons, in two indepen- 
dent one-dimensional bands 



H ° - E £ uk f lkJ,ka + E 



v k vk vk 



(4) 



with e v% = 



2tx b cos(fc^) + A — yti the spinon energy and 
2t\ cos (fey) + A the holon energy. Fixing x, the 
average number of holons per site (thus also the average 
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number of spinons) and using eq. (3), we can then self- 
consistently determine the values of A, \i, \^ an d X b - hi 
this way we incorporate effectively the spin-charge sepa- 
ration of the one-dimensional correlated electron system 
realized in the nearest-neighbor hopping Hamiltonian (1). 
We now add the next-nearest-neighbor intcrband hopping 
via two O 2p-orbitals as displayed in fig. lb). This weaker 
hopping term connects different one-dimensional electron 
systems and may be written in the slave-boson represen- 
tation as 



H' = 



a k.k' ,q 



(g b J)Kf . +h.c) (5) 

V q J lk+q<y lk 2k' J 2k'+qa J v ; 



with the structure factor gq = —4t' sm(q x ) sin(q y ) and t' 
the intcrband hopping integral. We now use the fact that 
between these systems only electrons can be transfcred, in- 
voking that a spinon and holon have to correlate to yield 
hopping. Thus we may interpret Ji! as an effective inter- 
action term introducing attractive coupling between the 
two subspecies with the tendency to recombine them to 
an electronic quasiparticle [13]. 

We now investigate the electronic spectrum using the 
retarded single-electron Green's function defined as 

G™'(q,t) = -e(t)({ % - CT (t), C t /f(j (0)}). (6) 

where the electron operators can be replaced by the spinon 
and holon operators, 



uqa 



-Y fK b 



(7) 



and Q(t) is the Heaviside step-function. In eq. (7), N is 
the number of lattice sites in one direction. 
Following standard methods we first calculate an auxiliary 
Green's function 

gf(q,t) = -eW^jW/^W,^ (0)}) (8) 

which can be evaluated in energy space using the equation 
of motion, 

Eg v /(q,E) = — y7iV # ft b \ 



H,bKf T _ 

vk vk-\-qa 



v qa 



(9) 



This equation of motion involves commutators of the form 




Fig. 2: Density log-plot of the spectral density function along 
the main symmetry axis for T = 0.005t showing the dispersion 
of the quasiparticles ('electrons'). The dotted lines mark the 
boundaries of the spinon-holon continuum. The actual quasi- 
particle band is thus strongly renormalized. The dashed, white 
horizontal line denotes the Fermi energy. 



The higher order Green's functions coming from (11) are 
treated within RPA to yield 

G(q, E) = G (q, E) + G (q, E)g(q)G(q, E) (12) 



where 



k-\-q,k 



is the bare Green's function resulting from the Hamilto- 
nian (4) without inter-band hopping and 



(15) 



Eventually, we find the RPA form of the single-electron 
Green's function in SBMFT 



G^(q,E) 



l-( gi yGl(q,E)Gl(q,E) 



(16) 



Results. — The total spectral density in momentum- 
space associated with the Green's function (16), 

S(q,E) = -Uv{lm{G™{q,E)]}, (17) 

7T 
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u ' vk J vk+qo\ k+q,k v k J vk+qa v ' 
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+ E 



is evaluated numerically where the hole concentration is 
x = 0.33 and the inter-band hopping is t'/t = 0.15 as 
parameters fitting to the de Haas-van Alphen Fermi sur- 
face [14]. For q along the x or the y axis in momentum 
space, gq vanishes and the denominator of (16) is equal to 1 
without any modifications to the original one-dimensional 
situation. In contrast, for q along the diagonal, \g^\ is 
maximal and the renormalization of the Green's function 
is strongest. Indeed the recombination of spinons and 
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Fig. 3: Comparison of the spectral density at q — (7r/2,7r/2) 
as a function of energy for the cases of low and high temper- 
atures with interband hopping and low temperature without 
interband hopping. 



holons into quasiparticles is here most pronounced, visible 
in the band hybridization and the evolution of quasiparti- 
cle weight as we will see below. 

First, we examine the spectrum along the main symme- 
try axis in the Brillouin zone (see fig. 2). At high tempera- 
tures, this spectrum is represented by the two-particle con- 
tinuum of a spinon and a holon, whose boundary is marked 
by the dotted lines in fig. 2. At low temperatures the spec- 
tral weight shifts into two quasiparticlc bands with a sharp 
momentum energy relation. These bands have the shape 
which is expected from a simple tight-binding model tak- 
ing nearest- and ncxt-ncarcst-ncighbor hybridization into 
account. Note that we ignore the lifetime effect due to 
the scattering among the resulting quasiparticles. Thus, 
these quasiparticle states remain well defined even away 
from the Fermi energy. The shift of spectral weight to the 
quasiparticlc bands is further illustrated in fig. 3 which 
shows a cut of the spectral density at q = (tt/2,tt/2) for 
different temperatures and interband couplings. We see 
in this plot that the peak around E = — 0.3£ without an 
interband coupling splits into two well-defined peaks of 
which the one at lower energies stays even for higher tem- 
peratures while the one at higher energies gets washed 
out. 

Second, we further analyze the q-dependence as well as 
the temperature dependence of the quasiparticle spectral 
weight. Figure 4 shows density log-plots of the spectral 
function at the Fermi energy for different temperatures. 
At lowest temperatures (T = O.OOli) the two-dimensional 
Fermi surface is well defined by quasiparticles (fig. 4a)). 
With increasing temperature, the quasiparticle weight de- 
creases and the Fermi surface becomes gradually faint. 
The retreat of the quasiparticle weight is clearly away 
from the [100] and [010] directions towards the [110] direc- 
tion (figs. 4.b-d)). A similar behavior of the quasiparticle 
weight was observed in ARPES experiments [3]. How- 
ever, only the 7-band was studied in the [110] direction 
and thus, a complete comparison of the results is beyond 




Fig. 4: Density log-plots of the spectral density in momentum 
space at the Fermi energy for different temperatures. For T — 
O.OOli, we see a Fermi surface as expected for a hybridized, 
two-dimensional system. When increasing the temperature the 
peaks of the spectral function get smeared and eventually, for 
T = O.lt we find a picture similar to that of a system of two 
decoupled one-dimensional bands. 



the scope of our approach. 

In our calculation, quasiparticles occur first along the diag- 
onal momentum direction and extend slowly over a wider 
range of the emerging Fermi surface to form in an inter- 
mediate stage pieces of arc-like Fermi surfaces. To illus- 
trate the evolution of quasiparticles at the Fermi surface, 
in fig. 5 we show the angle dependence of the quasipar- 
ticle weight of the pocket around the M -point for differ- 
ent temperatures. This clearly demonstrates the gradual 
build-up of quasi particle weight as the inset of fig. 5 shows. 
Within our calculation the Fermi arcs appear only clearly 
for the Fermi surface around the M-point (ct-Fermi sur- 
face). While for the Fermi surface centered around the 
T-point the rise of quasiparticle weight does not involve 
such a well defined arc feature. 

Note that unlike in other slave-boson discussions, here the 
holon condensation is not involved in the formation of co- 
herent quasiparticles. The formation of the quasiparticles 
in the spectrum is caused by the inter-orbital hybridization 
which is taken into account in the RPA-like form of the 
Green's function (16). The picture obtained in our slave- 
boson mean field approach should agree on a qualitative 
level with other descriptions of the dimensional crossover 
starting out from a Luttinger-liquid in one dimension. The 
slave-boson mean-field theory carries the advantage of for- 
mal simplicity. 

c-axis transport. — For the transport along the c- 
axis we assume that the interlayer tunneling approach is 
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Fig. 5: Quasi-particle weights as calculated from the spectral 
density function in arbitrary units around the M-point as a 
function of the angle to the x-axis for different temperatures. 
When decreasing the temperature a growing arc of quasipar- 
ticle weight appears around the diagonals. To emphasize this, 
the inset shows the endpoints of the arcs for different temper- 
atures. 



justified due to the extremely large anisotropy of resistiv- 
ity. We introduce the tunneling Hamiltonian 

H T = ^V V (T~ v 'fK .6,6+/ + h.c.) 

(18) 

where operators with momentum k belong to one plane 
and with p to the other. The matrix elements are 



^'=4t''^^,cos fc 



'sin^sin^ (19) 



for hopping from one band to the same and from one band 
to the other, respectively. This form results from the anal- 
ysis of the possible matrix elements between different lay- 
ers taking symmetry aspects into account. The total tun- 
neling current can be expressed as the sum of a coherent 
and an incoherent tunneling current 



hot — *coh + ^ir 



(20) 



Depending on the temperature range we expect one or 
the other part to dominate the transport. The maximum 
in c-axis resistivity in Sr2Ru04 is around T = 130A', 
for t f» 0.3cV [14] this corresponds to a temperature of 
T w 0.05*. We therefore have to study a low-temperature 
region with T <C 0.05* and a 'high'-temperature region 
with T > 0.05*. 

In the low temperature limit, we saw that the quasipar- 
ticle weight of the coherent spinon-holon pairs is growing 
with decreasing temperature and we expect the main con- 
tribution to the tunneling current [15] to come from the 
coherent component, 



vv' kp 



l S^k, 0)S»(p,0). (21) 



Here, Sl(u) denotes the spectral density of the band v 
in the lower (upper) layer. Using the spectral density as 
calculated from eq. (17) this indeed yields an increasing 
resistivity as is shown in fig. 6a), because of the dimin- 
ishing of quasiparticle weight upon growing temperature. 
Note that in this discussion we ignore any other source of 
temperature dependence. 

In the high-temperature limit, we see from figs. 4 and 
5 that coherent quasiparticles that could tunnel between 
layers have disappeared. Therefore, the quasiparticle pic- 
ture based on paired spinons and holons as used for the 
discussion of the tunneling current above is inappropriate. 
We have to take a different approach whereby spin and 
charge degrees of freedom are independent. In the tunnel- 
ing process of an electron, a spinon and a holon (charge 
carrier) have to be transfered simultaneously giving rise 
to an incoherent transport [16]. We start with the Fermi 
Golden Rule to estimate the transfer probability of holons, 



F 



2- 



J2\^o\n T m\ 2 S(E - E' + eV) (22) 



with cj)Q the many-body ground-state, <f>' any excited state 
and Eq and E' the corresponding energies. Decoupling 
the holons and the spinons, we find for the total rate of 
holon tunneling 



F 



^i 2 E « 

k,k' .q 



4\q)\l-nf(q)]x 



xn%\k')[l + nf{k)]S(u k - u k , + eV).(23) 

which requires that a spinon moves in the direction oppo- 
site to the holon. Here, np, B ^(k) denotes the Fermi (Bose) 
distribution with momentum k in the layer i, \Tj_\ 2 is the 
(momentum independent) tunneling matrix element. The 
factor of 8 in the above formula is composed of a factor 
of 2 for spins and a factor of 4 because of the different 
bands involved. Since the distributions are the same on 
both layers, we can write the net tunneling rate as 



2tt 



IT, 



F° — 

1 net 



deDp(e)nF{t)[l — np(e)] 
eV)[n B {u + eV) -n B (w)](24) 



where D F r B \ is the density of states of the spinons 
(holons). The integral over the Fermi distribution func- 
tions yields a factor of DpT with Dp the density of states 
of the spinons taken to be constant at the spinon Fermi 
level. Taking the density of states of the holons, D Bl to 
be constant too, the integral over the Bose distributions 
can be written as 

(D B ) 2 J dw{n B {u + eV) — n B (ui)} » D B n%K B eV, (25) 



with the compressibility of the hard-core holons, 
chemical potential, fj, and n B the holon density. 
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Fig. 6: The resistivity in the low- and high-temperature limit 
as calculated from the eqs. (21) and (27) .a) As expected, the 
decreasing quasiparticle weight in the low-temperature limit 
results in an increase in the resistivity as the temperature is 
increased, b) In the high-temperature limit, the tunneling pos- 
sibility of the spinons dominates the picture resulting in a de- 
creasing resistivity. 

Based on eq. (24) the resulting current can now be esti- 
mated as 

I = ^{2v)\T^\ 2 D F D B Tn 2 B K B {T)V (26) 
leading to the conductivity 

2 

o = 8j(2K)\T ± \ 2 D F D B Tn%K B (T)-^. (27) 

The compressibility of the holons, k b (T), corresponds to 
that of spinless Fermions and is constant for tempera- 
tures much smaller than the characteristic energy of holons 
which is of the order ■ Thus, the temperature depen- 
dence of a in the range of interest is dominated by T due 
to the spinon phase space contribution which is increasing 
with growing temperature. Therefore we find a decreasing 
resistivity in the temperature range above T ~ 0.05t. 
In addition to the c-axis resistivity, the analysis of the 
magnetic susceptibility displays signatures of a dimen- 
sional crossover, as magnetic excitations due to nest- 
ing features, pronounced in one dimension, are weakened 
when lowering the temperature. A behavior of such kind 
was experimentally reported in [17]. 

Conclusion. — In our discussion we have interpreted 
the regime change of the c-axis transport in Sr2Ru04 in 
terms of a dimensional crossover (from one to two dimen- 
sions) for the electronic states in the basal plane. We have 
shown that the slave-boson approximation can give a qual- 
itative understanding and provides interesting insight into 
the change of transport properties. In particular, our re- 
sults lead to a picture where with decreasing temperature 



a Fermi surface appears through gradually extending arcs 
of finite quasiparticle weight. This kind of Fermi surface 
evolution could be accessible to more specific ARPES ex- 
periments. In view of the progressive emergence of co- 
herent quasiparticles the difference between the metallic 
and insulating behavior can be understood as a change 
between dominant coherent to dominant incoherent tun- 
neling as temperature is raised. 
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